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The axes of gyroscopes experimentally define non-rotating frames. But 
what physical cause governs the time-evolution of gyroscope axes ? Starting 
from an unperturbed, spatially flat FRW cosmology, we consider cosmological 
vorticity perturbations (i.e. vector perturbations, rotational perturbations) at 
the linear level. We ask: Will cosmological rotational perturbations drag the 
axis of a gyroscope relative to the directions (geodesies) to galaxies beyond the 
rotational perturbation? We cast the laws of Gravitomagnetism into a form 
showing clearly the close correspondence with the laws of ordinary magnetism. 
Our results are: 

1) The dragging of a gyroscope axis by rotational perturbations beyond the 
H radius {H = Hubble constant) is exponentially suppressed. 

2) If the perturbation is a homogeneous rotation inside a radius significantly 
larger than the H radius, then the dragging of the gyroscope axis by the 
rotational perturbation is exact for any equation of state for cosmological 
matter. 

3) The time-evolution of a gyroscope axis exactly follows a specific average 
of the matter inside the H radius for any equation of state. 

In this precise sense Mach's Principle follows from cosmology with Einstein 
Gravity. 
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1 Introduction and Conclusions 



The spin axes of gyroscopes define our local nonrotating frame experimentally. But what 
physical cause governs the time-evolution of gyroscope axes? Newton invoked "absolute 
space" . Mach [Q] made the hypothesis that the nonrotating frame is determined (caused) 
"in some way" by the motions of all masses in the universe. He did not know any 
mechanism, any force, which would have this effect. 

We shall show that General Relativity, specifically Gravitomagnetism, in the context of 
the Friedmann- Robertson- Walker (FRW) cosmology, spatially flat, with rotational per- 
turbations (treated as linear perturbations) makes predictions in total agreement with 
Mach's hypothesis as formulated above. 

Thirring in 1918 analyzed the partial dragging of inertial frames inside a rotating 
infinitely thin spherical shell of uniform surface mass density and total mass M. He found 
that general relativity in the weak field approximation and in first order in the angular 
velocity i^sheii predicts an acceleration of test particles in the interior corresponding to a 
Coriolis force, which could be eliminated by going to a reference frame which is rotating 
with Vt' = /drag^sheii- the Weak field approximation, G^M <^ i?, he obtained the 
dragging fraction /drag = 1- Lense and Thirring made the corresponding 

analysis outside a rotating star. 

Einstein, after initially being inspired by Mach's ideas, concluded in 1949 |0: "Mach 
conjectures that inertia would have to depend upon the interaction of masses, precisely as 
was true for Newton's other forces, a conception which for a long time I considered as in 
principle the correct one. It presupposes implicitly, however, that the basic theory should 
be of the general type of Newton's mechanics : masses and their interactions as the original 
concepts. The attempt at such a solution does not fit into a consistent field theory, as will 
be immediately recognized." We shall show in section 4, how this apparent contradiction 
with Mach's conjecture is resolved in General Relativity, specifically in Gravitomagnetism. 

Goedel in 1949 @] presented a model which contradicts Mach's principle, but since it 
contains closed time-like curves, it is not relevant for physics. 

Brill and Cohen in 1966 again considered a very thin rotating spherical shell to 
lowest order in the rotation frequency, but now to all orders in the mass M of the shell, 
i.e. they treated geometry as a perturbation of Schwarzschild geometry. They found that 
space-time is fiat throughout the interior and that the interior Minkowski space is rotat- 
ing relative to the asymptotic Minkowski space. For large masses, whose Schwarzschild 
radius approaches the shell radius. Brill and Cohen found that the induced rotation inside 
approaches the rotation of the shell, i.e. the dragging of the inertial frame inside by the 
rotating matter of the shell becomes perfect. 

A cosmologically more relevant example is a rotating ball of dust and a gyroscope 
at the center. In a weak-field approximation with a Minkowski background this is the 
same as a superposition of a sequence of Thirring shells, which gives /drag of the order of 
GArMbaiiZ-Rbaii- The Friedman equation gives GnP of order of H"^ = R]f, if the curvature 
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term is not dominant. Hence /drag is of order of -Rbaii/-^?f) valid for i?baii Rh- Hence the 
dragging factor approaches the order of 1 only for i?baii Rh, but in this case the weak 
field approximation on a Minkowski background breaks down. See also the discussion of 
Mach's principle in Misner, Thorne, and Wheeler 0. 

C. Klein analyzed a thin spherical shell with empty interior embedded in a Fried- 
mann universe. The shell is required to follow the expanding motion of the surrounding 
cosmic dust. Rotations are treated to first order in the angular velocity. In his treat- 
ment the dragging coefficient tends to one, if nearly "the whole mass of the universe" is 
concentrated in the thin shell. 

At a conference on Mach's Principle in 1993, edited by Barbour and Pfister the 
question "Is general relativity with appropriate boundary conditions of closure of some 
kind perfectly Machian?" was put to the participants at the end of the conference. The 
vote was "No" with a clear majority (page 106 of f^). 

In this paper we discuss a realistic cosmological model and use cosmological per- 
turbation theory including super-horizon perturbations (instead of perturbations around 
Minkowski space or around the Schwarzschild solution), we use realistic cosmological mat- 
ter (instead of matter with a contrived energy- momentum-stress tensor), and we analyze 
the most general cosmological perturbations in the vorticity sector (instead of toy models 
like rigidly rotating thin shells). We consider a background FRW cosmology, spatially 
fiat, with linear vorticity perturbations, and we obtain the following specific results: 

1) The dragging of gyroscope axes by rotational perturbations beyond their H radius 
is exponentially suppressed, where H = Hubble constant. 

2) For a homogeneous rotation of cosmological matter inside a perturbation radius Rp, 
the dragging of the axis of a gyroscope at the center approaches exact dragging 
exponentially fast as Rp increases beyond the H radius. This holds for any equation 
of state for cosmological matter. 

3) A gyroscope axis exactly follows a specific average of the energy fiow in the universe 
with an exponential cutoff outside the H radius for any equation of state. 

In this precise sense we have shown that for a spatially fiat FRW Universe with ro- 
tational perturbations (added at the linear level) Mach's Principle on nonrotating frames 
follows from General Relativity (without the need to impose any boundary condition of 
closure). 

2 Vorticity Perturbations 

The vorticity a; of a velocity field v is defined by a; = curl v. Every velocity field can 
be decomposed into a potential fiow, curl v = 0, and a vorticity fiow, div v = 0. The 
simplest vorticity perturbation is a spherical region which rotates with a homogeneous 
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angular velocity around an axis through its center, v = f2e x r . A more realistic 
cosmological perturbation is one with a Gaussian cutoff, fl{r) = Qq exp {—r^/Rp). In the 
cosmological context a perturbation in a spherical region with a velocity field regular at 
the origin is much more appropriate than the vorticity field familiar from hydrodynamics 
with cylindrical symmetry and a vortex-line singularity in the center. 

The most general vorticity field v can be given by an expansion in vector spherical 
harmonics X^^m(6',0), ref. [@: 

Xe,m{OA) = LYe,m{eA), L = -tfxV , (1) 
vir) = 5]A,,^(r) X,,„(^,0) + 5]curl [5,,^(r) 0)] . (2) 

The simplest case is i = 1, which can be oriented such that m = 0: The A-solution is a 
homogeneous rotation at each value of r (see the simple examples above), which is called a 
"toroidal" field. The -B-solution is called a "poloidal field" . If the velocity field is toroidal, 
the vorticity field is poloidal and vice versa. These two configurations are familiar from 
ordinary magnetostatics, where curl B = An J. 



3 Cosmological Perturbation Theory 

For superhorizon perturbations we need a general relativistic treatment. This was devel- 
oped by J. Bardeen in his Phys. Rev. article of 1980 (ref. ^^). In linear perturbation 
theory one has three decoupled sectors of perturbations: scalar (density perturbations), 
vector (vorticity or rotational perturbations), tensor (gravitational wave perturbations) 
under rotations in 3-space. In the vector sector all perturbation quantities must be built 
from 3-vector fields with zero divergence, i.e. from pure vorticity fields. 
The perturbation 6goo is a scalar under 3-rotations, hence 

5900 = . (3) 

Therefore the slicing of space-time by space-like hypersurfaces of fixed time Ej has mea- 
sured times between slices for At = 1 (lapse) unperturbed. 

In contrast it is not useful (although possible) to insist on a time-orthogonal foliation 
(hence time-orthogonal coordinates). Rather we allow that the lines x = const are not 
orthogonal on the hypersurfaces of constant time Sf. It is useful to introduce the concept 
of a field of fiducial observers (FIDO's), which are located at fixed values of x, ref. [^. 
The shift-3- vector field (3 is defined as the 3- velocity dv/dt of the FIDO relative to the 
normals on S^. 

One can easily show ||10| that it is always possible to find a gauge transformation 
(coordinate transformation) such that the new coordinates give 

5(3) g.^ = . (4) 
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Therefore the 3-geometry of is unperturbed. We take the FRW background spatially 
flat, hence has the geometry of the Euclidean 3-space. In the gauge of eq. (^) has 
comoving Cartesian coordinates, 

ds^ = -dt^ + a^{t){dx'f + 2al3i dx' dt . (5) 

Geodesies on St are straight lines on our choice of chart. In all other gauges the coordinates 
of Euclidean 3-space are such that if the coordinates are Cartesian at one time, at all other 
times the coordinate lines (e.g. the x-axis) are no longer geodesies, they are getting wound 
up relative to geodesies. 

Asymptotically (beyond the region of the perturbation) we have FRW geometry and 
FRW coordinates, hence distant galaxies (galaxies beyond the rotational perturbation) 
are at fixed x (comoving coordinates). Therefore the coordinate axes in our gauge stay 
fixed relative to distant galaxies, i.e. the coordinate axes are geodesies on and nonro- 
tating relative to distant galaxies. In the context of rotating black holes such coordinates 
are called "star-fixed coordinates". For brevity of notation we shall employ this short 
term, although we always mean "coordinates fixed to galaxies beyond the rotational per- 
turbation" . 



4 The Laws of Gravitomagnetism for 
"star"-fixed FIDO's 

By what dynamical mechanism can cosmological matter in rotational motion far away 
influence the spin axes of gyroscopes here? By the laws of general relativity, more precisely 
gravitomagnetism. Since we work to first order in vorticity perturbations, we must analyze 
weak-field gravitomagnetism. We cast the laws of weak-field gravitomagnetism into a form 
which shows clearly the close correspondence with the familiar laws of ordinary magnetism 
(electromagnetism). We introduce and work with the gravitomagnetic field Bg and the 
gravitoelectric field Eg measured by fiducial observers (FIDO's) ref. |TI|, where the FIDO 



is a crucial concept for our formulation of the laws of gravitomagnetism. Note that for 
observers who are free-falling and non-rotating relative to gyroscopes, i.e. for inertial 
observers, there are no gravitational forces. Eg = 0, Bg = 0, by the equivalence principle. 

Our choice is to work with "star"-fixed FIDO's (more precisely FIDO's fixed by galax- 
ies beyond the rotational perturbation). "Star" -fixed FIDO's are defined to stay at fixed 
X in the "star"-fixed coordinate system, and (as in general) eQ(FIDO) = 'u(FIDO). Fur- 
thermore the spatial unit vectors ej of the FIDO's LONB (local orthonormal basis) point 
towards fixed distant galaxies. 

The 3-momentum measured by FIDO's is denoted by p = ej. Hats on indices refer 
to LONB's, the components jP are quantities directly measured by FIDO's. Note that 
this 3- vector lies in the hyperplane orthogonal to u(FIDO), not in the hyperplane S^. 
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The time- derivative dp/dt measured by FIDO's for a free- falling test particle gives the 
operational definition of the gravitational force, 



{dp/dt) bee fall =■ Fg ■ (6) 

— * 

The gravitoelectric field Eg is operationally defined via quasistatic test particles, 

4 := (^] ■ m 



m 



quasistatic particle 



Eg — g is the gravitational acceleration relative to the FIDO (i.e. measured by a FIDO) 
for a free-faUing particle. Compare with FIDO's on the surface of the earth (at fixed r), 

— * / 2 — * ~* 

who measure g = —9.8 m/s . The gravitomagnetic field Bg is defined via the first-order 
term in the velocity of the test particle (we put c = 1), 

— » — * 

{Fy)l^t order in = X Bg) . (8) 

After having given the operational definitions of Eg and Bg, which are valid for any 
choice of FIDO's, and having fixed our choice of FIDO's ("star" -fixed), we can derive the 
laws of gravitomagnetism. The results are: 

1) Equation of motion for test particles. We give it for vorticity perturbations of 
Minkowski space (for simplicity in this conference report), 

dp/dt = e (Eg + vx Bg) , (9) 

where e is the energy measured by the FIDO. This law is valid for free- falling 
particles of arbitrary velocities (e.g. photons), but only for "star"-fixed FIDO's. 
For other FIDO's there will be terms bihnear in the velocity. Our law should be 
compared to the Lorentz law of electromagnetism dp/dt — q{E + v x B). 

2) Equation of motion for the spin of a test particle or gyroscope at rest with respect 
to the FIDO. 

1 - - 

dS/dt = --BgXS. (10) 

This gives an angular velocity of precession of the spin axis relative to the axes of 
the FIDO (which in our case is relative to the direction to distant galaxies), 

^precession 2 ^9 ' (■^■^) 

This should be compared to the analogous precession of the intrinsic angular mo- 
mentum (spin) of classical charged particles in an electromagnetic field, i^precession = 
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+ 2 ^ B. Note again that for an inertial observer dp/dt = 0, dS/dt = 0. Every- 
thing depends on the choice of the field of FIDO's. Using eqs. (§) and (pH]) we 

obtain Fg — 2'm.(f2p]-ecession 

X v), i.e. the Coriohs force, where we must remember 
that fiprecession is minus the rotation velocity of the FIDO relative to the gyroscope 
axis. Since we are working to lowest order in Q the centrifugal term vanishes. A 
homogeneous gravitomagnetic field can be transformed away completely by going to 
a rigidly rotating coordinate system, i.e. physics on a merry-go-round in Minkowski 
space is equivalent to physics in a homogeneous gravitomagnetic field. This is analo- 
gous to the equivalence between physics in a homogeneous Newtonian gravitational 
field and physics in a linearly accelerated frame. 

3) Relation between Eg, Eg and the shift vector (3. The gravitational metric perturba- 
tions 5 g^y, i.e. the gravitational potentials for vorticity perturbations in "star"-fixed 
coordinates are given by the shift vector /3, which is connected to Eg and Eg by 

Eg = -dt0, Eg = VX^ . (12) 

We see that the shift vector (3 plays exactly the same role for gravitomagnetism 
as the electromagnetic vector potential A for electromagnetism. Therefore we can 
write (3 = Ag, and call it the gravitomagnetic vector potential. 

4) Einstein's Cqj equation for weak-field gravitomagnetism: 

A) In Minkowski background 

V xBg = -IQn Gjv Je , (13) 

where is the energy current, and in linear perturbation theory = {p + p)v 
with V = Compare Einstein's Gq^ equation to the Ampere-Maxwell equation 

in electromagnetism, V x E — dt E = +471 Jq, where Jq is the charge current. In 
weak-field gravitomagnetism for "star"-fixed FIDO's the Ggj equation is a constraint 
equation (i.e. no time derivatives, an equation at one moment of time with the 
same form as Ampere's law of magnetostatics), while the Ampere-Maxwell equation 
contains a time-derivative and is not a constraint equation. 

B) In a FRW background eq. (0) gets an extra term, 

V X Bg-4H (3 = -IGn Gn Je ■ (14) 

/^From H = —AttGnIp + p) see that H < for p > —p, therefore we define the H 
radius by i?^ = {—H)~^, and we define /i^ = —AH = AR~^. We insert the vector 

potential Ag = (3, we use div Ag = and curl curl Ag = —A Ag, hence we obtain 
by Fourier transformation 

(^Shys + /^') ^g = -16n GnX- (15) 
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We see that the new term on the left-hand side, {—AH (5) resp (/i^ Ag), dominates for 

— * 

superhorizon perturbations. The solution is the Yukawa potential for the source Jg, 

= /3(r,t) = -AGn I d'r'X{r',t) ^M-j^\r-r'\) ^ ^^^^ 
J \f — 

analogous to the formula for ordinary magnetostatics except for the exponential 
cutoff. 



5 Mach's Principle 



The Yukawa potential in eq. (|1^) has an exponential cutoff for |r — r'| > ^, hence we 



reach our first important conclusion: The contributions of vorticity perturbations beyond 
the H radius are exponentially suppressed. 

Our second result concerns the exact dragging of gyroscope axes by a homogeneous 
rotation of cosmological matter out to significantly beyond the H radius (for the expo- 
nential cutoff to be effective). This holds for any equation of state. This is easily seen 
from Einstein's Gqj equation in fc-space for superhorizon perturbations, /Cpj^yg ^ {—H), 

-4H(3 = -IQtiGnJs, Je = {p + p)va^id, H = -47iGm{p + p). (17) 

All the prefactors cancel, and we obtain /3(x) = —vauidix). With f2gyroscope = = 
— |(V X P) from eqs. ([TI|) and (|12]), and with Qnuid = ^ ^^fluid) we obtain 

— * — * 

^gyroscope ^fluid- (1^) 

This proves exact dragging of gyroscope axes here by a homogeneous rotation of cosmo- 
logical matter out to significantly beyond the H radius. 

Our third result concerns the most general vorticity perturbation in linear approxi- 
mation, and it states what specific average of energy flow in the universe determines the 
motion of gyroscope axes here at r = 0, 

^gyroscope('" 0; ^) 2 ^9^^ 0) ^) 

2Gn{p + p) fd^r\ \{l + pr) e"^"! [f x v (r, t)] . (19) 



The right-hand side of this equation is the gravitomagnetic moment of the energy current 
distribution on a shell [r, r + dr], analogous to the magnetic moment of an electric current 
distribution, integrated over r with the given weight function. This is the lowest term, 
the / = 1 term, in the multipole expansion for robs = and rgource > 0. Higher multipoles 
do not contribute to flgyroscope at r = 0. Only the toroidal velocity field for / = 1 of eq. ([l|) 
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contributes in eq. (JT^, i.e. a term corresponding to a rigid rotation velocity f2(r) for each 
shell (r, r + dr). 

Eqs. (|T6|) and (|I9[) are equations at fixed time, the precession of a gyroscope axis here 
is determined by the rotational motion of the masses in the universe at the same cosmic 
time. This shows that Einstein's objection to Mach's principle of 1949 |^ is not valid for 
weak gravitomagnetism. Furthermore, because of the exponential cutoff, which follows 
from Einstein's equations in the cosmological context, there is no need to impose "an 
appropriate boundary condition of some kind" . 

Is exact dragging a measurable effect in principle? There seem to be two problems: 1) 
To obtain exact dragging (within some given experimental error) the observational cutoff 
radius must be significantly larger than the H radius (for the exponential cutoff to be 
effective). We have defined f^gyroscopo and fifluid relative to distant galaxies, i.e. galaxies 
beyond the rotational perturbation. But we cannot "see" galaxies beyond the H radius. 
2) The Gqj equation for "star"-fixed FIDO's is an equation at fixed time, our results hold 
at fixed time, but an experimental test today with telescopes here on earth cannot see the 
galaxies at the same cosmic time, it sees galaxies on the past light cone. The solution to 
these problems is an experimental test in principle: Take data locally by many different 
observers all over x space, also beyond the H radius. Collect and patch together the data 
after a few Hubble times. This can be done as a computer experiment, where the resulting 
test only involves our measurable quantities. This computer experiment is not necessary, 
since we have given the proof (at the level of linear perturbation theory). 

We conclude that in a spatially flat FRW Universe with rotational perturbations 
(treated at the linear level) Mach's Principle on nonrotating frames follows from Gen- 
eral Relativity. 
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